Anisotropy and Ising-like transition of the S — 5/2 two-dimensional Heisenberg 

antiferromagnet Mn-formate di-Urea 
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Recently reported measurements of specific heat on the compound Mn-formate di-Urea (Mn-f- 
2U) by Takeda et al. [Phys. Rev. B 63, 024425 (2001)] are considered. As a model to describe 
the overall thermodynamic behavior of such compound, the easy-axis two-dimensional Heisenberg 
antiferromagnet is proposed and studied by means of the pure quantum self-consistent harmonic 
approximation (PQSCHA). In particular it is shown that, when the temperature decreases, the 
compound exhibits a crossover from 2D-Heisenberg to 2D-Ising behavior, followed by a 2D-Ising- 
like phase transition, whose location allows to get a reliable estimate of the easy-axis anisotropy 
driving the transition itself. Below the critical temperature T N = 3.77 K, the specific heat is well 
described by the two-dimensional easy-axis model down to a temperature T* = 1.47 K where a 
T 3 -law sets in, possibly marking a low-temperature crossover of magnetic fluctuations from two to 
three dimensions. 

PACS numbers: 75.10.Jm, 05.30.-d, 75.40.-s, 75.40.Cx 



The quantum Heisenberg antiferromagnet (QHAF) on 
the square lattice is one of the most widely investigated 
magnetic models because of both its fundamental the- 
oretical properties and the existence of many real com- 
pounds whose magnetic behavior is ruled by a spin-spin 
interaction properly described by the 2D QHAF Hamil- 
tonian. Due to their layered structures, such compounds 
have an intralayer exchange integral J which is much 
larger than the interlayer one, J'. This ensures a 2D 
thermodynamic behavior to persist down to a certain 
crossover temperature which could be naively estimated 
of the order of J'S 2 . However, most of the above men- 
tioned compounds display a phase transition towards an 
ordered phase at a critical temperature of the order of 
JS 2 . 

If one assumes a fully isotropic intralayer interaction, 
which cannot induce any finite temperature transition 
by itself, this experimental finding may be explained by 
noticing that strong correlations between spins on a sin- 
gle layer act as an effective amplifier of the interlayer in- 
teraction, in that the coupling between neighboring spins 
on different layers drags in a similar coupling the sur- 
rounding spins within a distance of the order of the 2D 
magnetic correlation length*. However, most of the ex- 
perimentally observed phase transitions do also display 
features which suggest a persistent 2D behavior, and are 
not compatible with the above sketched mechanism. On 
the other hand, such experimental observations may be 
explained by allowing the intralayer interaction to be 
anisotropic; the existence of some anisotropic coupling 
may hence be invoked, and a more detailed analysis de- 
veloped. 



In this paper we consider the quasi-two-dimensional 
real compound Mn(HCOO) 2 -2(NH 2 ) 2 CO (Mn-formate 
di-Urea, or Mn-f-2U), which is is a remarkable realiza- 
tion of a S — 5/2 2D QHAF in a wide temperature 
region. From specific-heat and susceptibility measure- 
ments 2 ^* the intralayer spin-spin coupling J turns out 
to be J = 0.68 ± 0.04 K, while J' is orders of magnitude 
lower. Characteristic features of a phase-transition, such 
as a sharp peak in the specific heat, are also observed at 
T N = 3.77 ± 0.02 K and suggested^ the existence of an 
easy-axis anisotropy in the intralayer spin-spin coupling. 
By proton NMR measurements it has been established in 
Ref. |5j that the interlayer coupling is much smaller than 
the anisotropy, which has been estimated to be of the 
order of 10 -2 J. 

We hence propose, as a model for the thermodynamics 
of Mn-f-2U, the S = 5/2 QHAF on the square lattice 
with easy-axis exchange anisotropy (EA-QHAF), whose 
Hamiltonian reads 
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where i = (i\, i 2 ) runs over the sites of a square lattice, d 
connects each site to its four nearest neighbors, J > is 
the antiferromagnetic exchange integral and /i 6 [0, 1) is 
the easy-axis anisotropy parameter. The spin operators 
S? (a = x, y, z) are such that \S\ 2 = S(S + 1) and obey 

Classical^ and quantum* Monte Carlo simulations pre- 
dict such model to display a 2D-Ising phase transition 
at a finite temperature T^S*, [i], which continuously de- 
creases as [i increases, finally vanishing for fi —> 1, i.e., in 
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the isotropic model. The picture in the case of Mn-f-2U is 
hence that its 2D-Ising transition at T t immediately trig- 
gers the observed*^* phase transition to 3D long-range 
order: therefore, we assume Tj = T N . 

We use the pure-quantum self- consistent harmonic ap- 
proximation (PQSCHA)&, a semiclassical method which 
reduces the expressions of quantum statistical averages 
to effective classical-like ones, where temperature- and 
spin-dependent renormalization parameters appear; the 
thermodynamics of the effective model can then be stud- 
ied by means of classical techniques, like classical Monte 
Carlo simulations. Besides a uniform additive term that 
does not affect the evaluation of statistical averages, the 
PQSCHA effective Hamiltonian for the EA-QHAF reads* 
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where S = S + 1/2 and Sj = (s?,s|,sf) are classical 
unit vectors. J e g(t, S, n) < J and fj, e g (t, S, fi) > fi are 
the renormalized exchange and anisotropy of the effective 
classical model, respectively. Hereafter the dimensionless 
temperature t — T/JS 2 is used. 

This method was applied to the EA-QHAF, and, in 
particular, to the study of the fi- and S'-dependence of 
the Ising critical temperature, which is shown in the 
phase diagram t Y (S, fx) reported in Ref. U Using that 
phase diagram, the anisotropy value corresponding to 
the experimental critical temperature of Mn-f-2U, t 1 = 
0.616±0.039, is n = 0-981±g;g^|, the large error being due 
to the experimental uncertainty on J. This anisotropy 
estimate is close to that of the compound Rb 2 MnF4, a 
S = 5/2 EA-QHAF we have extensively investigated in 
Ref. I lCt namely /i = 0.9942. Moreover, performing a new 
set of classical MC simulations with a slightly different 
value of /i, more fit to the case of Mn-f-2U, is, at present, 
too large an effort compared to the expected benefits. 
As a matter of fact, due to the poor knowledge of J, any 
choice of the value of /i in the interval estimated above is 
affected by a high degree of arbitrariness. We will hence 
hereafter use the already known PQSCHA curves** for 
/i = 0.9942; owing to the proximity between the two 
anisotropics (i.e., /i = 0.981+QQ2g estimated for Mn-f-2U 
from our phase diagram*, and \i = 0.9942 corresponding 
to our previous theoretical data—), we expect to be able 
to single out the characteristic features of the experimen- 
tal data by comparison with the theoretical curves. 

Focusing our attention on recent measurements* of the 
magnetic specific heat of Mn-f-2U, we make a twofold 
comparison of such data with our predictions for the EA- 
QHAF. The first comparison, shown in Fig.[I] is made be- 
tween the measured specific heat of Mn-f-2U and the the- 
oretical one with /j, = 0.9942 as a function of the dimen- 
sionless temperature. Besides these data, two PQSCHA 
curves for the S = 5/2 isotropic model are also reported: 
the dashed one is obtained within the same approxima- 
tion level as for the anisotropic model, while the dotted 
one is obtained by an improved version of the PQSCHA 
scheme^* that was only applied to the isotropic case; the 
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FIG. 1: Specific heat vs t = T/JS 2 , for S=5/2. Mn-f-2U 
experiments (squares)*, EA-QHAF with /i = 0.9942 (solid 
line)*, isotropic QHAF (dashed line), isotropic QHAF with 
improved PQSCHA (dotted line)—. Error bars are due to the 
experimental uncertainty on J for Mn-f-2U. Inset: zoom of 
the Heisenberg-Ising crossover region, with few representative 
error bars. 



latter was already compared to the experimental data by 
Takeda et alA, so it is included in order to avoid mistakes. 
Both the theoretical specific heat for the EA-QHAF with 
H = 0.9942 and the measured specific heat of Mn-f-2U are 
seen to behave as that of the isotropic QHAF for temper- 
atures well above t 1 . The agreement with the isotropic 
behavior, however, gets worse once the anisotropy starts 
to drive the anisotropic systems towards the Ising critical 
regime. The temperature where the anisotropic (solid) 
curve deviates from the isotropic (dashed) one to form 
the Ising spike marks a crossover between the Heisenberg 
and the Ising regime. This crossover, the position of the 
peak, and the whole behavior below the critical temper- 
ature depend on the anisotropy value /i. We stress that 
if we had calculated the PQSCHA curve for the value 
fi = 0.981 estimated from the PQSCHA phase diagram 9 , 
the position of the theoretical peak would coincide with 
that of the experimental one. 

In order to emphasize the common features of the the- 
oretical and the experimental curves in the Ising regime 
we scale their abscissa t with the actual critical temper- 
atures, i.e., t l = TJJS 2 = 0.616 for the experimental 
data, and t : = 0.575 for the theoretical ones, making the 
peaks match as shown in Fig. [5] In this way, remark- 
able agreement between experiment and theory shows 
up not only for the shape of the transition peak but also 
for the whole temperature region between t/t 1 = 1, and 
t*/tj ~ 0.42, where the experimental data for the com- 
pound Mn-f-2U depart from the PQSCHA curve. 

These results show that the compound Mn-f-2U, be- 
sides behaving as a 2D EA-QHAF above and in the criti- 
cal region, displays such behavior in its specific heath also 
below i : , down to a much lower temperature t* ~ 0.24. 
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FIG. 2: Specific heat vs t/t,, with t r = TJJS 2 = 0.616 for 
the experimental data, and t : = 0.575 for the theoretical line. 
The dashed line is a best fit of the first 8 experimental data 
points to a t 3 law; other symbols are as in Fig. 

This suggests that below t l , even though the weak inter- 
layer coupling drives the system into a 3D ordered phase, 
magnetic fluctuations, upon which the specific heat di- 
rectly depends, remain confined within the layers, be- 
ing those between different layers incoherent. The 2D 
picture of the ordered phase eventually breaks down at 
t = t* . Below this point the experimental data for Mn- 
f-2U are much better described by a simple i 3 -law; such 
law characterizes the low-t region of a 3D quantum an- 
tiferromagnet, with the main contribution arising from 
linear excitations. The temperature t* could thus be in- 
terpreted as marking the dimensional crossover from 2D 
to 3D behavior, corresponding to the onset of coherence 
in the magnetic excitations propagating perpendicularly 
to the layers. However, at this stage, the microscopic 



mechanism responsible for the observed low-t behavior 
cannot be firmly identified; nonetheless, the above pic- 
ture seems to be the simplest one that accounts for the 
experimentally observed behavior. 

In conclusion, we have compared the theoretical find- 
ings based on the PQSCHA for the 2D EA-QHAF with 
recent specific heat data on Mn-f-2U. The comparison re- 
veals the existence of a crossover from a high-temperature 
2D-Heisenberg regime to a critical 2D-Ising regime that 
triggers the phase transition at T N = 3.77 K observed 
in Ref. @. In the ordered phase, the specific heath still 
behaves as that of a 2D EA-QHAF down to the tem- 
perature T* = JS 2 t* ~ 1.47 K, where 3D behavior sets 
in. The success of the EA-QHAF model in describing 
the specific heat of Mn-f-2U in a wide temperature range 
makes trustworthy our prediction for the value of the 
anisotropy parameter, fi — 0.981, as obtained above by 
the sole knowledge of t l = T l / (JS 2 ). The above value of 

implies the anisotropy to be about 1.9% of the exchange 
energy. Accounting for the large uncertainty upon the 
estimated J, however, the actual anisotropy could range 
from 0.6% to 4.5%. Anyhow, the above estimate of the 
anisotropy strength is compatible with the one made by 
Kubo et al. in Ref. 0. Experiments on Mn-f-2U aimed at 
a more precise determination of both J and /i, as well as 
at the characterization of magnetic excitations, such as 
inelastic neutron scattering ones, come now to be highly 
desirable. 
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